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ABSTRACT 



DUDLEY KNOX UBf?',Ry 

naval postgraduate S-HOn, 
MONTEREY, CAL, PORN, A Co 



A software package was developed to provide the inputs for a 
passive single- sensor acoustic tracker. This required a computation 
of sound propagation velocity and resolving high- resolution frequency 
information from raw ?*coustic data. A selected frequency was then 
automatically tracked from an up-doppler to a down-doppler condition. 
This frequency information coupled with associated bearing information 
was then used to provide continuous target fixes. 
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L INTRODUCTION 



Continuous fixing of an underwater sound source by passive means 
has traditionally been carried out by using multiple sensors. Single 
sensor fixing requires varying power coefficients, frequencies, and/or 
bearing changes and usually produces an estimated position with an 
accuracy tliat varies greatly with the type of sensor and processing 
utilized and the skill of the evaluator. An accurate single sensor track 
would be advantageous in that positional errors could be minimized with 
only one buoy or array required, providing greater flexibility in that a 
series of sensors are not required to be monitored, and a decrease in 
buoy expenditures could be expected. To provide the necessary frequency 
data, a software package was developed that computed the sound propaga- 
tion velocity, detected a discrete frequency signal and automatically 
tracked that signal from a full up-doppler through a full down-doppler 
situation. This frequency output combined with associated bearing in- 
formation, and a value for the sound propagation speed with an estimate 
of the speed of the target provides svifficient input data for single sensor 
fixing. 
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II. THE TRACKING PROBLEM 



A. STATEMENT OF THE PROBLEM 

Chapter 3 of [1] devclopes the required input quantities for passive, 
single-sensor fixing. These values are an estimate of target speed, the 
velocity of sound propagation in the medium, the doppler shifted fre- 
quency and a bearing from the sensor to the target. The estimate of 
target speed v/ould have to be determined by tactical considerations at 
the time and location of signal acquisition, i. e. , target transiting or on- 
station, signal strength, detection ranges, etc. The bearing would be 

provided by a standard DIFAR buoy or an array. In this problem the 

o 

standard deviation of bearing measurement noise was varied froJii 5 lu 
30^. Methods of providing the sound propagation velocity and the doppler 
shifted frequency were then devised. 

B. DEVELOPING THE INPUTS 

Equation 1.1. of [2] provides a good approxiination of sound speed 

in sea water as a function of temperature, salinity and pressure. It is 

repeated here as equation 1. 

, 2 

c(T, S, z) = 1449+4. 62T-. 054T +1. 3(S-35)-i. 017z (1) 

where c is in meters per second, T is in degrees centigrade, S is in 
parts per thousand, and z is in meters. Converting this equation so 
that c is in yards per second, T is in degrees Farenheit and depth is in 
feet, yields equation 2. 



11 



( 2 ) 



c(T, S, z) = 1476. 1+3. 97T-. 018T^ + 1. 42(S-35)+. 0057z 
If the target depth is known, the median depth of the target and hydro- 
phone should be used. Otherwise, use the hydrophone depth. Salinity 
would normally be set to 35ppt unless operating in water of a known 
different salinity. The resultant output from this equation would be set 
equal to the velocity of sound propagation, VP, and used as an input 
value to the target tracker program. 

The frequency requirements called for a high- resolution value 
which dictated a relatively long time window. A twenty second window 
was used which resulted in a . 05Hz FFT resolution. This resolution 
was further enhanced by interpolation using the relative values of the 
power coefficients associated with the discrete frequency being tracked. 
In order to maintain contact on the same discrete frequency line, a 
frequency tracker algor itlim was devised. This provided the necessary 
frequency data for the target tracker program. 
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III. SOFTWARE IMPLEMENTATION 



A. SOUND SPEED 

This process was a simple matter of implementing the equation for 
sound speed in seawater, equation 2, Providing the chosen depth, ob- 
served temperature, and salinity results in a best estimate of the cictual 
sound speed in the operating area. Updating this value would be necessary 
as the operating area or environmental conditions changed. 

B. FREQUENCY 

1. Parameters for Frequency Comp utations 

In order to provide realistic frequency inputs, a magnetic tape 
was obtained that contained analog, raw acoustic signals representing 
ambient noise plus an underwater sound emitter with various discrete 
frequencies. The doppler shift of a discrete frequency signal is directly 
proportional to the rest frequency of the signal. Thus doppler shift 
measurements arc more easily obtained from higher frequencies. 

However, signal attenuation in the medium is much higher at higher 
frequencies. At present, frequencies in the 150Hz to 300Kz range will 
yield usable values v/ith realizable frequency resolution measurements. 

The highest signal frequency of interest determines the required sam- 
pling rate of the data. The sampling rate and the time window length 
determine the required dimension of the data input matrix. With these 
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considerations in mind, the following parameters were chosen in digitiz- 
ing the data: 

'Sampling Rate of 512 samples /second yields Nyquist frequency 
of 256Hz; 

Input data low-pass filtered to 250Hz to prevent aliasing errors; 

Time record length of 20 seconds for FFT resolution of . 05 
Hz/Bin; 

Hanning windov/ function applied in time domain to reduce 
sidelobes associated with rectangular data windows. 

2. Power Spectrum 

With the above chosen data parameters, the input matri>: to the 
POWER Subroutine''" [S] was dimensioned to 20,480 bins. This cor- 
responds to two bins per data point (1 cosine term and 1 sine term, for 
10,240 data points) or 10,240 Complex Fourier coefficients. The real 
data values were loaded into the odd (cosine or real) bins and zero’s 
were provided for the even (sine or imaginary) bins. On return from 
the FFT Subroutine [6] to the POWER Subroutine, the cosine term is 
squared and added to its associated sine term squared thus representing 
the power at some discrete frequency. This discrete frequency is an 
integer multiple of the FFT processing resolution (the reeiproca] of the 
time record length). The first bin, Y(l), contains the zero frequency 
power coefficient, the second bin contains the . 05 Hz power coefficient, 
and so on in .05 Hz steps. Thus bins Y(l) through Y(5001) now contained 
the power coefficients from zero to 25 0 Hz. . 



Subroutine STATS and POWER were written by Arfman, J. F. , Jr. , 
and were adapted for use in this program. 
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3. Search Windov/ and Background Noise Suppression 



V/'ith over five thousand power coefficients available for track- 
ing, it is convenient to limit the range of search to those that offer the 
best data for doppler tracking. A signal having favorable characteristics 
is first identified, i. e. , frequency above 100 to 150 Hz, apparently 
stable on a one Hz resolution processor, and of sufficient strength to 
maintain a favorable signal-to-nois e ratio for a majority of the time. 

Then a search window is designated by the adjustable parameters FL 
and SRBW, where FL is the lower frequency limit of the search v/indow 
and SRBW is the search window bandwidth. Both parameters must be 
some integer multiple of the FFT resolution and would usually be kept 
at whole numbers for convenience. The total number of power coefficients 
in the search window would equal N=20(SRBW) and are designated YA(1) to 
YA(N). 

Ideally when the frequency to be tracked (FTRK) has been iden- 
tified, FL and SRBW should be chosen such that FTRK is equal to FL 
plus (1/2)SRBW and SRBW is as large as is reasonable to keep the ratio 

of signal power in the window to the total power in the window much less 
L M 

than one. If S YA(K) « S YA(J) for L signal bins and M noise 
K=1 J=1 

bins with L“idv4=N, then the summation as I goes from one to N of YA(I)/N 
is a reasonable approximation of the ambient noise level eiround FTRK [4], 
[ 5 ]. This calculation is performed and set equal to YNOIS. Then YKOIS, 
or some function of it dependent on tlie signal-to-noi sc ratio, is sub- 
tracted from each power coefficient in the search window. Under the 
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assumption that the noise power spectrum is constant (white noise) 
within SRBV/, only the signal power spectrum remains in the window. 

Frequency Tracking to Follov/ the Doppler Shift 

The frequency selected to be tracked will be the strongest sig- 
nal in the search window when the first transform is computed. The 
bin number of that signal is saved as TBIN(JJ). In each succeeding 
transform the strongest signal (largest coefficient) in the search window 
is designated TBIN(JJ) and is compared to TBIN(JJ-l). A bin shift 
greater than 10 is determined to be a new or different discrete signal. 
The new power coefficients around the last center bin arc then examined 
for discrete signal energy. The maximum value found would be con- 
sidered the signal energy being tracked if its bin number is equal to the 
last center bin number plus or minus 1, and its magnitude is greater 
than a chosen threshold minimum. In this problem this value v/as set 
to three times the ambient noise value. If the maximum coefficient 
found is more than one bin from the last center bin, the old center bin is 
examined. If its coefficient is larger than at least one of its sidelobe 
bin values and at least half as large as the ma.ximum found and larger 
than the chosen threshold minimum value, the old signal bin is con- 
sidered to still contain the signal being tracked. If any of these test 
fail, the new maximum found is compared to the chosen minimum 
value. If the new maximum is larger, it is accepted as the signal being 
tracked; if not the contact is considered lost. Figure 1 possibly 
clarifies the above process. 
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LOAD NEW POWER COEFFICIENTS 
OF BINS AROUND TBIN(JJ-l) AND 
FIND NEW MAXIMUM COEFFICIENT. 
RESET TBIN(JJ) TO NEW VALUE AND 
COMPUTE NEW BNSHF. 







THMIN = A CI-IOSEN THRESHOLD MINIMUM VALUE 
FIGURE 1. FREQUENCY TRACKER DECISION FLOW. 
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Once the bin containing the largest portion of the signal energy 
being tracked has been established, the frequency of that signal is 
refined by computing the first moment of the three bin power ’’hump*' 
around it. If the detected signal is an integer multiple of the J'FT 
resolution, adjacent bin spill over will be minimum and the refining 
calculation would make no correction to the frequency of the bin 
associated with the maximum power. When the detected signal is not 
an integer multiple of the FFT resolution the refining calculation adjusts 
the frequency to a mid- resolution value corresponding to the first mo- 
ment of the hump. See Figure 2. The frequency corresponding to bin 
number K is stored in the YL matrix as YL(K). The frequency refining 
equation used was 

FTRK{JJ) = YL(K)'fRESOL^^(YA(K-IT)-YA(K-l))/(YA(K+l) 
+YA(K)+YA(K-1)) . (3) 

where 

FTRK=the discrete frequency being tracked, 

RESOL^the resolution of the FFT processing, 

YA(M)=the power coefficient of the Mtli bin. 

Letting YL(K)=220. 25Hz and using the sample values from figure 2, 

/ 1-4 

FTRK(JJ) = 220. 250Hz + . 05Hz ‘ ( 

V 1+7+4 

= 220. 250Hz + . 05IIz ‘ (-.25) 

= 220. 250Hz - .0125Hz 
= 220.2375Hz 

220.238Hz (rounded) 




18 



YA(K) 

(POWER) 



FIGURE 2. 







4 




FTRK=220. 238Hz 

/ i 

I !<“=» YL(K)=220.250Hz 

' \ 

/ I 
I 1 
/ \ 

/ I 
\ 

\ 

V 

V 
\ 




. K-1 K+1 

K 



FREQUENCY BINS 



RESOLVING MID-E.ESOLUTION FREQUENCIES WHEN THE 
DETECTED SIGNAL IS NOT AN INTEGER MULTIPLE OF 
THE FFT RKSOLUTION. 
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Without this calculation the returned value of the frequency being tracked 
would have been 220. 250Hz. A five bin computation was found to offer 
no significant improvement over the three bin computation due to the 
increased probability of noise or adjacent discrete signal influence. 

C. DESIGNATING PROBLEM VALUES 

In computing the sound propagation speed for the following examples, 
nominal values were used because the true values were not available. 
Values used v/ere: 

Temperature 48^F 
Depth = 60 ft 

Salinity = 35 ppt 

The sound emitter on the data tape appeared to make a slow turn 
into the buoy for the first 4 minutes of the tape and then maintain a 
constant course at a speed of approxiinately 4.5 knots. The closest- 
point-of-approach (CPA) to the hydrophone was approxima.tely 300 to 700 
yards and 20 minutes of data was utilized with CPA at about the 11 to 12 
minute point on the tape. A 20 Hz search window bandwidth was used 
with several differejit lower frequency limits on the window. It was 
found that one transform per minute coiTifortably contained the tracked 
frequencies within the 21 bin tracking v/indow in all cases for a slow- 
moving, near-field emitter. 
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D. DATA OUTPUT 



1. Run 1 (See Table 1) 

V/ith FL set to 40 Hz and SRBW at ZO Hz, the rest frequency 
of the tracked signal was approximately 50. 13 Hz, The total doppler 
shift was only 0. 15 Hz and was too low to provide good tracking data. 

2. Run 2 (See Table 2) 

For this run FL was set to 220 Hz with SRBW still at 20 Hz. 

Two strong discrete frequency signals were within the search windov/ 
bandwidth, one’s rest frequency at 233. 87 Hz and the other at 223.22Hz. 
The maximum power coefficient changed betv/een these two signals 
several times. In the initial transform the tracker locked on to the 
lower of the two frequencies and a continuous track was maintained even 
though the signal faded out and then came back. This run gave the fre- 
quency tracker a good check-out in that (a) the maximum coefficient in 
the search window had considerable variations between discrete signals 
(b) a noise spike occurred in the tracker window at time=430 (FILE^^ZZO) 
and (c) the tracked signal faded completely out at time=1090 (FILE=550) 
and was regained one minute later. The total doppler shift was approx- 
imately 0. 6 Hz and provided good data for doppler tracking. The computer 
output for this run is included in this report in its entirety. 

3. Run 3 (See Table 3) 

For this run FL was set to 225 Hz and SRBW remained at 20 Hz. 
The only significant signal in the search window was the 233. 87 Hz signal 
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Time 


Center 


Ratio: 


Frequency- 


(seconds) 


Bin 


Signal/Noise 


Tracked 


10. 0 


204 


155/66 


50. 155 


70. 0 


205 


1759/38 


50. 190 


130. 0 


205 


1619/59 


50. 185 


190. 0 


205 


2001/54 


50. 203 


250. 0 


205 


1620/39 


50. 180 


310. 0 


205 


1527/49 


50. 192 


370. 0 


204 


878/40 


50. 197 


430. 0 


205 


247/42 


50. 151 


490. 0 


205 


9320/78 


50. 178 


550. 0 


205 


13280/68 


50. 191 


610. 0 


204 


56746/254 


50. 168 


670. 0 


203 


38878/144 


50. 105 


730. 0 


203 


42170/187 


50. 085 


790. 0 


203 


47640/227 


50. 081 


850. 0 


202 


6941/55 


50. 060 


910. 0 


202 


9654/71 


50. 063 


970. 0 


202 


18841/125 


50. 053 


1030. 0 


202 


17558/111 


50. 057 


1090. 0 


202 


17546/104 


50. 054 


1150. 0 


202 


14377/87 


50. 050 



Rest Frequency approximately 50. 13 Hz 



TABLE I. 

FREQUENCY TRACKER OUTPUT WITH FL= 50 Hz, SRBW=20 Hz. 
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Time 


Center 


Ratio: 


Frequency 


(seconds) 


Bin 


Signal/Nois e 


Tracked 


10. 0 


70 


69/6 


223. 437 


70. 0 


71 


112/5 


223. 481 


130. 0 


71 


61/5 


223. 495 


190. 0 


72 


52/5 


223. 556 


250. 0 


71 


254/6 


223.489 


310. 0 


70 


218/6 


223. 473 


370. 0 


69 


490/9 


223. 420 


430. 0 


69 


64/8 


223. 407 


490. 0 


70 


743/24 


223. 434 


550. 0 


71 


2131/41 


223. 488 


610. 0 


69 


10429/68 


223. 396 


670. 0 


62 


4795/50 


223. 063 


730. 0 


60 


2554/26 


222. 960 


790. 0 


61 


2119/16 


222. 995 


850. 0 


60 


322/14 


222, 938 


910 . 0 


59 


245/10 


222. 882 


970. 0 


59 


60/4 


222, 898 


1030. 0 


59 


112/6 


222. 886 


1090. 0 


60 


12/5 


Lost Track 


1150. 0 


59 


115/5 


222. 901 



Rest Frequency approximately 223.22 Hz 



TABLE II. 

FREQUENCY TRACKER OUTPUT WITH FL=220 Hz, SRBW=20 Hz. 
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y. 






■■ 











Time 


Center 


Ratio: 


(seconds) 


Bin 


Signal/Nois e 


10. 0 


183 


61/4 


70, 0 


184 


42/4 


130. 0 


185 


21/4 


190, 0 


185 


77/5 


250. 0 


184 


167/5 


310, 0 


184 


160/5 


370. 0 


183 


829/7 


430. 0 


184 


708/8 


490. 0 


184 


3295/19 


550. 0 


184 


4824/32 


610. 0 


183 


4669/30 


670. 0 


176 


6049/29 


730. 0 


174 


906/13 


790. 0 


174 


766/7 


850. 0 


172 


1082/11 


910. 0 


172 


1020/8 


970. 0 


172 


27/4 


1030. 0 


172 


357/5 


1090. 0 


172 


135/5 


1150. 0 


171 ■ 


23/5 



Rest Frequency approximately 233, 87 Hz 



TABLE III, 

FREQUENCY TRACKER OUTPUT WITH FL=225 Hz, 
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Frequency 
T racked 



234. 100 
234. 151 
234. 193 
234. 224 
234. 169 
234. 163 

234. 102 
234, 144 
234. 132 
234. 130 
234. 079 
233. 756 
233. 644 
233. 634 
233. 563 
233.558 

235. 563 
233. 545 
233, 530 
233. 525 



SRBW=20 Hz. 



and this was tracked with only a minimal amount of decision require- 
ments around CPA, The total doppler shift was approximately 0. 65 Hz 
and provided good data for doppler tracking. 

E. TARGET TRACKING 

Since no bearing information was available corresponding to the 
frequency data, the IBM 360 Computer was used to generate simulated 
real time bearings. This was accomplished by programming a pseudo 
*'true track” into the computer so that the target speed and the CPA time 
were matched as v/ell as possible. The initial starting position was 
X=500 yards, Y = -1050 yards at time T(l) = 190 seconds. The heading 
and speed used were 90 degrees and 4. 5 knots where zero degrees 
corresponds to the positive X axis. The computer thexi computed ”true 
bearings” to the target and added random noise values of a programmed 
standard deviation, SA, to these values. Thus the true track plotted 
by the computer is only a best estimate of the targets relative position 
from the buoy. 

The X-Y Filter developed by Mitschang [l] was used to process the 
data and develop the target track. To initialize the filter and commence 
tracking, the measured bearing must have changed by at least three 
standard deviations of the measurement noise and the doppler shifted 
frequency must have decreased from its value at the start of the initializa- 
tion process. A modification to the filter was attempted so that while 
initializing, an increasing frequency would be detected as a inancuvering 
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target. Time, bearing, and frequency values occurring with the maxi- 
mum frequency detected would then be retained as the initial values to 
be used when the bearing shift and non-increasing frequency criteria 
were again satisfied. This would delay the start of the tracking process 
but would prevent data from a maneuvering target being used to initialize 
the state equations based on a constant heading, constant speed target. 
However, indexing problems were encountered and were not resolved 
in time to include the modification in this report. 

The following difference equations developed in' Chapter 3 of [1] 
were used to compute the initial range, heading, and rest frequency of 
the target. 



R = -(VP* DELF- DELT)/(f(last)' (DEL0) 



2 



(4) 



0si=0ave-18O - arcsin((- VP- DEEF)/(f(last)- DEL0- Vsi)) 



(5) 



F0i=Fave [l - (vsi/VP)cos(0si - 0ave)] 



( 6 ) 



where 



VP =speed of sound propagation 



DELF =doppler shift of the frequency 



BELT =time required to meet initialization criteria 



f(last) =last frequency measured for initialization 



DEL0 ^change in bearing from sensor to target 



0si =^initial estimate of target heading 



0ave ^average bearing for the initialization period 



Vsi ^initial estimate of target speed 
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FOi =initial estimate of the rest frequency 

Fave =average frequency for the initialization period 
Then based on a constant heading, constant speed target, the initial 
filter state equations are: 



X(l) 


=X position = R cos Gave 


(7) 


X(2) 


=X velocity = Vsi cos Gsi 


(8) 


X(3) 


= Y position = R sin Gave 


(9) 


X(4) 


= Y velocity = Vsi sin Gsi 


(10) 


X(5) 


= Rest frequency = FOi 


(11) 



The initial covariance matrix values are computed by the ^‘direct 
method** developed in Chapter 4 of [l] and an extended Kalman filtering 
technique was used to refine the estimates of the target’s position, 
heading, speed, and rest frequency as the problem progresses. 

Since the target was maneuvering for the first four minutes of the 
tape, initialization of the target tracker was commenced at time 
T(4) = 190 seconds, with the values from Table 3 as input data. The 
average initialization time, AVKJ, varied v/ith the standard deviation of 
the bearing measurement noise, SA. This time varied from about 4 
minutes for SA=5 degrees to about 8 minutes for SA=30 degrees. 
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IV. RESULTS OBTAINED FROM REAL FREQUENCY INPUTS 



Onee the frequency trackei’ has locked on to a stable but doppler 
shifted frequency, the inputs required of the operator are an estimate of 
the rest frequency, FO, and an estimate of the target*s speed, AVSK. 
These.are ^lpriming’‘ values used in computing the initial heading of the 
target and affect the ultimate accuracy of the computed track. A third 
variable affecting the initialization is the standard deviation of the bear- 
ing error, SA. As SA increases, the filter takes longer to start tracking 
and has more data on which to base its original estimates. 

The program was run for various combinations of the three variables 
FO, AVSK, and SA, The true value for FO was approximately 233. 87 IIz 
and for AVSK was 4, 5 knots. An estimated FO greater than the detected 
frequency causes a ”down-doppler solution for the initial course. Con- 
versely an estimate of FO less than the detected frequency results in an 
‘*up-doppler ” solution for the initial course. Each combination of var- 
iables was looped so that one hundred runs were made with that set of 
inputs and the averaged filtered positions were plotted. The symbol 
on each plot represents the computed position of the target at time T(l) 
plus one-half of the initialization period. The range is computed by 
Equation (4) and the bearing is the average bearing for the initialization 
period. The true points from whieh the noisy bearijig measurements 



28 



were made are plotted with the symbol Averaged filter positions 

are plotted with an "X". The ascending order of precedence on each 
plot is +, X. Figures 3 through Z4 show the performance of the system 
with real, high resolution frequencies as inputs. F(DET) signifies 
initial frequency detected. 
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FIGURE 3. 
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FIGURE 7. 
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FIGURE 8 
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FIGURE 9. 
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AVSK> TRUE AVSK, SA = 15°. 
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FIGURE 13. FILTER OUTPUT FOR EST. FO F(DET), EST. 
AVSK> TRUE AVSK, SA = 15°. 
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FIGURE 14. FILTER OUTPUT FOR EST. FO >- F(DET), EST. 
AVSK> TRUE AVSK, SA = 15°. 
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FIGURE 17. FILTER OUTPUT FOR EST. FO< F(DET), 
EST. AVSK> TRUE AVSK, SA = 10°. 
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FIGURE 18. FILTER OUTPUT FOR EST. FO > I'(DET), 
EST. AVSK > TRUE AVSK, SA = 10°. 
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FIGURE 19. FILTER OUTPUT FOR EST. FO F(DET), 
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FIGURE 20. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THANF(DET), EST. AVSK < TRUE AVSK, AND 
SA = 5°. 
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FIGURE 22. FILTER OUTPUT FOR EST. FO < F(DET), EST. 
AVSK> TRUE AVSK, SA = 5°. 
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V. CONCLUSIONS 



The automatic frequency tracker developed in this report has been 
shown to be adequate for tracking a doppler-shifting frequency in a 
favorable signal- to-noise ratio environment. The width of the tracking 
•’window^wouklrLave .to be .increased for a higher speed target unless the 
transforms were computed at closer time intervals, A loss of signal 
would result in the last center bin frequency being retained as the signal 
being tracked unless a spurious signal with more power than the chosen 
minimum appeared in the window. Therefore, under marginal signal- 
to-noise ratio conditions, an operator would be required to monitor the 
frequency tracker output to insure that after a 'lost track" condition 
occurred, any "regain track" operation was indeed occurring on the 
previous signal being tracked. 

The procedure of using Equation (3) to extend the frequency resolu- 
tion of the output above that of the FFT processing is dependent upon 
having a reasonable approximation of the ambient noise level. Cortipvit- 
ing a new estimate in each transform results in higher power levels 
than the true ambient noise level being used when the target is near CPA. 
Since the ratio of two numbers, A/B, is not equal to the ratio (A-K)/ 
(B-K) where K is a constant, the subtraction of the corrupted estimate 
of ambient noise level causes an error in the valvie of the resolved 
frequency. However, in comparing the amplitudes of the signal and 



52 



noise estimates around CPA (from Tables I, II, and III), the signal is 
150 to 200 times as large as the corrupted ambient noise estimate. The 
resultant error on the resolved frequency is on the order of magnitude of 
10"^Hz and was considered negligible. As the range at CPA increases, 
the signal power level decreases and its corrupting influence on the 
ambient noise level approximation is decreased. Thus as long as the 
increased average power level in the window is caused by the target 
being tracked, no appreciable error occurs and the resolution of the 
overall system has been extended beyond the resolution of the FFT 
processing. 

With real frequency inputs and reasonable initial estimates of the 
rest frequency and target speed, the Mitschang Filter provided target 
tracks in excellent agreement with the known emitter locations and 
probable emitter locations. Thus a true single sensor fix appears 
possible using existing hardware and techniques. Further testing is 
needed where exact target locations are known and real measured bear- 
ings are available corresponding to the real frequency measurements. 
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OfIPUTER OUTPUT 
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T2DF =(T2 - TH2{KA))/RA0 

TCP = i »/TA\'(T10F) - 1./TAN(T2DF) 

9A3 A?X!KA) = 02/ ( S I ^>! ( T 1 DF ) *TDF ) 

TF(AB5tARX(l) - AR HALF ) , L T . AB S ( AP X ( 2 ) - ARHALF) ) ARHALF - ARX(l) 
IF( AE5 (ARX(2) - ARHALF ), LT , ABS ( ARX { 1 ) - ARHALF)) ARHALF = ARX{2) 
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